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Before the launch of the Fermi satellite only two classes of AGN were known to generate relativistic jets and 
thus emit up to the 7-ray energy range: blazars and radio galaxies, both hosted in giant elliptical galaxies. The 
first two years of observations by the Large Area Telescope (LAT) on board Fermi confirmed that these two are 
the most numerous classes of identified sources in the extragalactic 7-ray sky, but the discovery of variable 7-ray 
emission from 5 radio-loud Narrow-Line Seyfert 1 galaxies (NLSyls) revealed the presence of a possible emerging 
third class of AGN with relativistic jets. Considering also that NLSyls are typically hosted in spiral galaxy, 
this finding poses intriguing questions about the nature of these objects, the onset of production of relativistic 
jets, and the cosmological evolution of radio-loud AGN. Here, we report on a preliminary investigation of the 
properties of this sample of radio-loud NLSyl at MeV-GeV photon energies, utilizing the four-year accumulation 
of Fermi LAT data. In addition we briefly discuss some radio-to-gamma-rays properties of the 7-ray emitting 
NLSyl in the context of the blazar scenario. 



1. Introduction 



Only a small percentage of Active Galactic Nuclei 
(AGNs) are radio-loud, and this characteristic is com- 
monly ascribed to the presence of a relativistic jet, 
roughly perpendicular to the accretion disk. Accre- 
tion of gas onto the super-massive black hole (SMBH) 
is thought to power these collimated jets, even if the 
nature of the coupling between the accretion disk and 
the jet is still among the outstan ding open qu estions 
in high-energy astrophysics [e.g. iMeierl [2003 1. Cer- 
tainly relativistic jets are the most extreme expres- 
sion of the power than can be generated by a SMBH 
in the center of an AGN. These objects have a total 



the Unification model for AGNs and the evolution of 
radio-loud AGNs. 

Even if the physics necessary to explain the genera- 
tion, collim ation and evolut i on of the jet is still to be 
known [e.g. lBlandfordl[200ll l2008j , there is increasing 
evidence that the properties of jets in AGNs could be 
related to the properties of the acc retion flow that 
feeds the c entral SMBH. Recently, iGhisellini et al.l 
20081 12011 | suggested a transition between BL Lac 



bolometric luminosity of up to 10 



49-50 



erg 



lAckermann et al 

20l4 iBonnoli et al. 2011 



[e.g. 

with a 

large fraction of the power emitted in 7 rays. Be- 
fore the launch of the Fermi satellite only two classes 
of AGNs were known to generate these structures and 
thus to emit up to the 7-ray energy range: blazars and 
radio galaxies, both hos ted in giant elliptical galaxies 
Blandford k Reedll978j . The first 2 years of observa- 



objects and Flat Spectrum Radio Quasars (FSRQs) 
that can be mainly justified by the different accretion 
regimes: sub-Eddington in the first class of objects, 
leading to radiatively inefficient accretion flows and 
relatively weak jets, whereas near-Eddington in the 
second class, giving rise to bright disks and powerful 
jets. This scenario seems to be also in agreement with 
the separation of their "non-beamed" counterparts, 
Fanaroff-Riley I and Fanaroff-Riley II radio galaxies. 
In this context, the case of the peculiar radio-loud 
NLSyls has received increasing attention. 

NLSyl is a c lass of AGN identified by 



tion by Fermi-L AT confirmed that the ext ragalactic 7- 
ray s ky is dominated by these two classes [Nolan et al.l 
|2012| . However, the discovery by the Large Area 
Telescope (LAT) on board the Fermi satellite of vari- 
able 7-ray emission from a few radio-loud Narrow- 
Line Seyfert Is (NLSyls) revealed the presence of 
a possible third cla ss of AGNs with relativistic jets 
Abdo et al.l [2009a, b.cj. On the contrary, no radio- 



TUll and characterized by 
narrow permitted lines 
km s -1 ) emitted from the 
< 3, and a Fell bump 



quiet Seyfert galaxies were d etected in 7 rays until 
now [Ackermann et ajj[2012bj . This finding poses in- 
triguing questions about the knowledge of the devel- 
opment of relativistic jets, the disk/jet connection, 



lOsterbrock k Poggd 
their optical properties: 
(FWHM (H/3) < 2000 
broad line region, [ OIII] /H /3 

[for a review see e.g. lPoggell2000j . They also exhibit 
strong X-ray variability, steep X-ray spectra, substan- 
tial soft X-ray excess and relatively high luminosity. 
These characteristics point to systems with smaller 
masses of the central BH (10 6 -10 8 M Q ) and higher 
accretion rates (close to or above the Eddington limit) 
with respect to blazars and radio galaxies. NLSyl are 
generally ra dio-quiet, with only a small fraction of 
them [< 7%, iKomossa et alj|20"ol| classified as radio- 



eConf C121028 



4 th Fermi Symposium : Monterey, CA : 28 Oct- 2 Nov 2012 



loud, and objects with high values of radio-loudness 
(R > 100) are even more sparse (~2.5%), while 
~15% of quasars are radio-loud. In the past, several 
authors inferred from studies of non-simultaneous 
radio to X-ray data dif ferent possib i lities on the 
nature of these objects. iKomossa et al.l |2006| argued 
that radio-loud NLSyls could be some young stage of 
quasars, while I Yuan et~aT, [2008] claimed that some 
NLSyl might be high-freque ncy-peaked flat s pectrum 
radio quasars conjectured bv lPadovanil [2007] and not 
predicted by the blazar sequence, but no conclusive 
evidence has been prese nted. However, as show n with 
preliminary results in iFoschini et al.l [2009j |. some 
radio-loud NLSyls show similarities with FSRQs, 
while others are like BL Lacs. Moreover, as in the 
case of blazars and radio galaxies, there should be 
a "parent population" with the jet viewed at large 
angles. Th e first source of th is type could be PKS 
0558-504 iGliozzi et ail l2010j. Other candidates 



have been detected bv lDoi et ail 2012 ], Therefore, it 
is possible that NLSyls are a set of low-mass systems 
parallel to blazars and radio galaxies. 

The firm confirmation of the existence of rclativis- 
tic jets also in Seyfert galaxies, provided by the LAT 
detection of NLSyls, opened a large and unexplored 
research space for important discoveries for our knowl- 
edge of the AGNs, but brought with itself new chal- 
lenging questions. What arc the differences between 
this class of 7-ray AGNs and blazars and radio galax- 
ies? What is the origin of the radio-loudness? What 
are the parameters determining the jet formation and 
how it is possible to have the formation of a jet in 
a population of AGN mostly radio-quiet such as the 
NLSyls? Is there a limiting BH mass above which 
objects are preferentially radio-loud? How do these 
objects fit into the blazar sequence? Four years after 
the announcement of the detection of th e first 7-ray 
NLSyl by Fermi-L AT, PMN J0948+0022 [Abdo et all 
l2009aj . only a few indications about the nature of 
these objects have been obtained. In Section 2 we 
report the LAT data analysis and results of the 5 7- 
ray NLSyls over 4 years of Fermi observations, fo- 
cusing on the two flaring sources SBS 0846+513 and 
PMN J0948+0022 in Section 3. In Section 4 and 5 
we discuss the host galaxies and the jet formation, 
respectively, for the NLSyls. In Section 6 we report 
about searching new 7-ray NLSyls, while concluding 
remarks are presented in Section 7. 



2. The 7-ray view of NLSyl 

Up to now 5 radio-loud NLSyl galaxies have been 
detected at high significance by Fermi-L AT: PMN 
J0948+0022, PKS 1502+036, 1H 0323+342, PKS 
2004-447, and SBS 0846+513. Here we analyze the 
first four years of 7-ray observations of these sources. 
The LAT data reported in this paper were collected 



from 2008 August 4 (MJD 54682) to 2012 August 
4 (MJD 56143). During this time the LAT instru- 
ment operated almost entirely in survey mode. The 
analysis was performed with the ScienceTools soft- 
ware package version v9r27pl. The LAT data were 
extracted within a 10° Region of Interest (Rol) cen- 
tred at the location of the 5 NLSyls. Only events be- 
longing to the "Source" class were used. In addition, 
a cut on the zenith angle (< 100°) was also applied 
to reduce contamination from the Earth limb 7 rays, 
which are produced by cosmic rays interacting with 
the upper atmosphere. The spectral analysis was per- 
formed with the instrument response functions (IRFs) 
P7S0URCE_V6 using an unbinned maximum likelihood 
method implemented in the Science tool gtlike. A 
Galactic diffuse emission model and isotropic compo- 
nent, which is the sum of an extragalactic and in- 
strumental background were used to model the back- 
ground 1 . The normalizations of both components in 
the background model were allowed to vary freely dur- 
ing the spectral fitting. 

We evaluated the significance of the 7-ray sig- 
nal from the sources by means of the Test Statis- 
tics TS = 2Alog(likeliho od) between mode ls with 
and without the source [Mattox et al.l Il996j ] . The 
source model used in gtlike includes all the point 
sources from the 2FGL catalogue that fall within 
20° from the target source. The spectra of these 
sources were parametrized by power-law functions, 
dN/dE oc (E I E )~ r , where T is the photon index, or 
log-parabola, dN/dE oc (E/Eo)-"- 13 ^s(.b/e )^ wnere 
E is a reference energy, a the spectral slope at the 
energy Eq, and the parameter /3 measures the curva- 
ture around the peak. A first maximum likelihood was 
performed to remove from the model the sources hav- 
ing TS < 10 and/or the predicted number of counts 
based on the fitted model N prec [ < 3. A second maxi- 
mum likelihood was performed on the updated source 
model. The fitting procedure has been performed with 
the sources within 10° from the target source included 
with the normalization factors and the photon indices 
left as free parameters. For the sources located be- 
tween 10° and 20° from our target we kept the nor- 
malization and the photon index fixed to the values 
of the 2FGL catalog. We used a power-law model 
for PKS 1502+036 and PKS 2004-447, and a log- 
parabola for 1H 0323+342 and PMN J0948+0022, as 
in the 2FGL catalog. For SBS 0846+513 we use d a 
power-law model as in lD'Ammando et al.l [2012a] for 
the third year of Fermi observation. 

The results of the LAT analysis over the entire pe- 
riod are summarized in Table |U The average observed 
isotropic luminosity of the 5 objects in the 0.1-100 



1 http: / / fermi . gsfc. nasa.gov / ssc / data / access/lat /Background 
Models.html 
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GeV energy range spans between 10 44 erg s _1 and 
10 47 erg s _1 , a range of values typical of blazars. 
This could be an indication of a similar viewing an- 
gle with respect to the jet axis and beaming factor for 
the 7-ray emission between blazars and 7-ray NLSyls. 
In particular, SBS 0846+513 and PMN J0948+0022 
showed 7-ray flaring activ ity combined with a mod- 
erate spectral evol ution [D' Ammando et all l2012al 
iFoschini et al.|[20TT |. a behaviour already observed in 
bright FSRQs and low-synchr otron-peaked (L SP) BL 
Lacs detected by Fermi LAT [Abdo et a,l.ll2010j . The 
average photon index of the 5 NLSyls spans a range 
(r = 2.2-2-6) narrower than those observed in Mis- 
aligned AGNs [F = 1.9-2.7. iGrandi fc Torres! l2012j 
and similar to the average values obs erved for FSRQs 
and LSP BL Lacs [Nolan et al.ll2012 |. 

Figures [1] to [5] show the 7-ray light curves of the 5 
NLSyls for the period 2008 August 4 -2012 August 4 
using 3-month time bins. For each time bin the spec- 
tral parameters of the target source and all sources 
within 10° from it were frozen to the value resulting 
from the likelihood analysis over the entire period. If 
TS < 10 the value of the fluxes were replaced by the 2- 
<j upper limits. The systematic uncertainty in the flux 
is energy dependent: it amounts to 10% at 100 MeV, 
decreasing to 5% at 56 MeV, and increasing a gain 
to 10% above 10 GeV [Ackermann et all l2012aj . A 
different level of activity has been clearly observed in 
these five objects, with significant variability in SBS 
0846+513 and PMN J0948+0022. 



2009 



2010 



20 1 



2012 



3. Flaring NLSyls: PMN J0948+0022 and 
SBS 0846+513 

One of the key question is the maximum power 
released by the jets of radio- loud NLSyl, and for 
this reason 7-ray flaring activities from these sources 
are catalyzing a growing interest in the astrophysi- 
cal community. The first answers arrived in 2010 
July when PMN J0948+0022 underwent a 7-ray flar- 
ing activit y with a daily peak flux of ^lxlO -6 ph 
cm~ 2 s _1 IFoschini et alj|2011 |. The first spectral en- 



ergy distributions (SEDs) collected for the 4 7-ray 
NLSyls detected in the first year of Fermi opera- 
tion showed clear similarities with blazars: a double- 
humped shape with disk component in UV, physical 
parameters blaza r-like and jet powe r in the average 
range of blazars Ab do et al . 2009c]. The compari- 
son of the SED of PMN J0948+0022 during the 2010 
July flaring activity with that of 3C 273, a typical 
FSRQ, shows a more extreme Compton dominance in 
the NLSyl. The disagreement of the two SEDs can 
be due to the difference in BH masses and Doppler 
factor of the two jets. Further 7-ray flaring episodes 
from PMN J0948+0022 have been observed in 2011 



June iD'Ammando fc Ciprinil koilH and 2013 January 
[D'Ammando fc Orientill2013j . 
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Figure 1: Integrated flux light curve of PMN 
J0948+0022 obtained in the 0.1-100 GeV energy range 
during 2008 August 4 - 2012 August 4 with 3-month 
time bins. Arrows refer to 2-a upper limits on the source 
flux. Upper limits are computed when TS < 10. 



A strong 7-ray flare was observed also from SBS 
0846+513 in 2011 June, reaching an isotropic 7-ray 
luminosity (0.1-300 GeV) of ~10 4 8 erg s" 1 , compara- 



ble to that of the bright FSRQs [D'Ammando et al.l 
2012a]. Variability and spectral properties of SBS 
0846+513 in radio and 7 rays indicate a blazar-like 
behaviour. In addition, from the model-fitting of 4- 
epoch MOJAVE data of this source in 2010-2012 we 
found that the jet components are separating with an 
apparent velocity of (10.9±1.4)c. This value suggest 
the presence of boosting e ffect as well as in blazars 
[D'Ammando et alj|2012b) . The power released dur- 
ing the flaring activity and the apparent superluminal 
motion velocity are strong indications of the presence 
of a relativistic jets as powerful as those in blazars, 
despite the lower BH masses. 



4. Host galaxy 

The discovery of a relativistic jet in a class of AGN 
usually hosted in spiral galaxies such as the NLSyls 
was a great surprise. Unfortunately only very sparse 
observations of the host galaxy of radio-loud NLSyls 
are av ailable up to now a nd t he sample of object s stud- 
ied by iDeo et all |2006| and IZhou et al.l [2006| had z 
< 0.03 and z < 0.1, respectively, including both radio- 
quiet and radio-loud sources. 

Among the 7-ray emitting NLSyls detected by 
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Table I Gamma-ray characteristics of radio-loud NLSyl galaxies detected by Fermi-LAT. 



Source Redshift Flux (E > 100 MeV) Photon index Curvature TS L 7 

z xlO -8 ph cm" 2 s" 1 T/a xlO^ergs" 1 



1H 0323+342 0.061 3.52±0.42 

SBS 0846+513 0.5835 2.63+0.22 

PMN J0948+0022 0.585 11.32+0.43 

PKS 1502+036 0.409 4.14+0.39 

PKS 2004-447 0.24 1.37+0.33 



2.62+0.13 0.63+0.17 199 0.01 

2.18+0.05 - 658 2.95 

2.37+0.05 0.23+0.03 2287 10.02 

2.60+0.06 - 309 1.25 

2.56+0.14 - 49 0.11 
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Figure 2: Integrated flux light curve of SBS 0846+513 
obtained in the 0.1-100 GeV energy range during 2008 
August 4 - 2012 August 4 with 3-month time bins. 
Arrows refer to 2-cr upper limits on the source flux. 
Upper limits are computed when TS < 10. 
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Figure 3: Integrated flux light curve of PKS 1502+036 
obtained in the 0.1-100 GeV energy range during 2008 
August 4 - 2012 August 4 with 3-month time bins. 
Arrows refer to 2-cr upper limits on the source flux. 
Upper limits are computed when TS < 10. 



Fermi LAT only for 1H 0323+342 Hubble Space Tele- 
scope (HST) and Nordic Optical Telescope observa- 
tions are available. These observations revealed a one- 
armed galaxy morphology or a circumnuclear ring, 
suggest ing two possibilitie s: the spiral arms of the host 
galaxy IZhou et aL 20071 or the residual of a galaxy 



Anton et al.l 20081 ] . On the other hand, no 



merging 

significant resolved structures hav e been observed in - 
stead by HST for SBS 0846+513 (Maoz et alJfl99^ . 
and no high-resolution observations are available for 
the remaining 7-ray NLSyls. Thus the possibility that 
the development of relativistic jets in these objects 
could be due to strong merger activity, not so usual 
in disk/spiral galaxies, cannot be ruled out. Further 
high-resolution observations of the host galaxy of 7- 
ray NLSyls will be fundamental to obtain important 



insights into the onset of production of relativistic jets. 



5. Radio-loudness and jet formation 

The mechanism at work for producing a relativistic 
jet is not still clear. In particular the physical parame- 
ters that drive the jet formation is still under debate. 
One fundamental parameter could be the BH mass, 
with only large masses t hat can allow effici ent forma- 
tion of a relativistic jet. ISikora et al.l [20071 ] suggested 
that AGN with M B h > 10 8 M have a radio loud- 
ness 3 orders of magnitude greater than the AGN with 
M B h < 3 x 10 7 M Q . According to the "modified spin 
paradigm" proposed, another fundamental parameter 
for the efficiency of a relativistic jet production should 
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Figure 4: Integrated flux light curve of 1H 0323+342 
obtained in the 0.1-100 GeV energy range during 2008 
August 4 - 2012 August 4 with 3-month time bins. 
Arrows refer to 2-cr upper limits on the source flux. 
Upper limits are computed when TS < 10. 
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Figure 5: Integrated flux light curve of PKS 2004-447 
obtained in the 0.1-100 GeV energy range during 2008 
August 4 - 2012 August 4 with 3-month time bins. 
Arrows refer to 2-cr upper limits on the source flux. 
Upper limits are computed when TS < 10. 



be the BH spin, with SMBHs in elliptical galaxies 
having on average much larger spins than SMBHs in 
spiral galaxies. This is due to the fact that the spi- 
ral galaxies are characterized by multiple accretion 
events with random angular momentum orientation 
and small increments of mass, while elliptical galax- 
ies underwent at least one major merger with large 
matter accretion triggering an efficient spin-up of the 
SMBHs. The accretion rate (thus the mass) and the 
spin of the BH seem to be related to the host galaxy, 
leading to the hypothesis that rel ativistic jets can de- 
velop only in elliptica l galaxy [e.g. lBottcher fc Dermerl 
120021 iMa rschcr 2009] . In this context the large radio- 
loudness of SBS 0846+513 could challenge this idea if 
the BH ma ss (8 2 xl0 6 -5.2xl0 7 M odot ) estimated by 
IZhou et al.1 [2005] is confirmed. 

We noted that BH masses of radio-loud NLSyls 
are generally larger with respect to th e entire sample 
of NLSyls [M BH «(2- 10) xlO 7 M^; iKomossa et all 
l2006l lYuan et all 12008} . even if still small when com- 
pared to radio-loud quasars. The larger BH masses of 
radio-loud NLSyls could be related to prolonged ac- 
cretion episodes that can spin-up the BHs. The small 
fraction of radio-loud NLSyls with respect to radio- 
loud quasars could be an indication that not in all 
of the formers the high-accretion regime laste d long 
enough to spin- up the central BH [Sikora 2009] . 



6. New 7-ray emitting NLSyl? 



In addition to the 5 NLSyls already observed 
new 7-ray emitting NLSyl could be detected ac- 
cumulating more and more Fermi-LAT data. For 
this reason we started from a list of 39 NLSyls 



repor ted in iKomossa et al.l I2006I . lOshlack et al.l 
I2001I. Iwhalen et all [2006} . HYuan et all [2008} . 
IZhou fc Wand [2002j with a radio-loudness R > 20 
and analyzed 4 years of LAT data to search systemat- 
ically for GeV emission from them. No significant (TS 
> 25) detection of new 7-ray NLSyls over the entire 
4-year period was obtained. 

Anyway the variability could be a key ingredient 
to take into account for a high-significance detection 
of these sources, as clearly showed by SBS 0846+513 
that after two years of very low 7-ray activity [and for 
this reason not included in t he Fi r st or Second Fermi- 
LAT catalog, lAbdo et al.l l2010L iNolan et all 12012} 
started an increase of the flux u p to the flaring ac- 
tivity detected in 2011 June- July [D'Ammando et all 
12012 a} . For this reason the next step will be an investi- 
gation of the activity of these candidates 7-ray NLSyl 
on different time scales (e.g. 1 month, 3 months, 1 
year). 
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7. Concluding remarks 

The presence of a relativistic jet in some radio-loud 
NLSyl galaxies, first suggested by their variable radio 
emission and flat radio spectra, is now confirmed by 
the Fermi-LAT detection of five NLSyls in 7 rays. 
The flaring episodes observed in 7 rays from SBS 
0846+513 and PMN J0948+0022 are strong indica- 
tions of the presence of a relativistic jet as power- 
ful as those of blazars. These two sources showed all 
the characteristics of the blazar phenomenon and they 
could be a relatively low mass (and possibly younger) 
version of blazars. Further multifrequency observa- 
tions of these and the other 7-ray emitting NLSyls 
will be fundamental for investigating in detail their 
characteristics over the entire electromagnetic spec- 
trum. The impact of the properties of the central en- 
gines in radio-loud NLSyls, which seem quite different 
of those of quasar and manifest in their peculiar opti- 
cal characteristics, on the 7-ray emission mechanisms 
is currently under debate. In addition, the detection 
of relativistic jets in a class of AGN usually hosted in 
spiral galaxies is very intriguing, challenging the the- 
oretical scenario of relativistic jet formation proposed 
up to now. The detection of new NLSyls in 7 rays by 
Fermi-LAT will be important for extending the sam- 
ple and better characterizing this new class of 7-ray 
emitting AGN. 
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